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1. Introduction 

Lattice QCD simulations enable us to calculate QCD thermodynamic quantities, such as the 
equation of state (EOS) in quark gluon plasma (QGP), non-perturbatively. The lattice results serve 
indispensable elements in understanding the nature of QGP in heavy Ion collision experiments, e.g. 
as inputs to the hydrodynamical space-time evolution of QGP. 

In the last decade, importance of dynamical quarks in EOS became gradually apparent. How- 
ever study of EOS with dynamical quarks requires one of the most expensive calculations in lattice 
QCD. It is sometimes hard to control lattice artifacts in the EOS calculations even with recent high 
performance computers. 

First of all, the lattice should be fine enough to keep lattice artifacts small. In the conventional 
fixed N t approach (N t is the temporal lattice extent), however, the lattice spacing a is large for 
low temperatures. Such coarse lattices not only cause lattice artifacts in the observables but also 
introduce uncertainties in the determination of the line of constant physics (LCP). Because the 
integral method [[I]] requires an integration from a low T, the uncertainties in the low T region 
affect EOS in the whole range of T. Therefore, it is important to keep the lattice fine in the whole 
range of T. On the other hand, in the high T region, small values of a there lead to small lattice 
volumes which may cause finite volume effects. 

In this study, we adopt the fixed-scale approach, in which we vary T by changing N t at a 
fixed a, by fixing all coupling parameters [0 0]. In this approach the problems mentioned above 
are in part resolved. Among others, we do not need to determine LCP at all: Since the coupling 
parameters are common to all temperatures, the condition to follow a LCP is automatically satisfied. 
Furthermore, when we borrow zero-temperature configurations from a large scale spectrum study, 
the lattice spacings are in the scaling region and are smaller than those used in conventional fixed 
N t studies around the transition temperature T c . In the high T region, on the other hand, while 
the lattice volume is kept fixed, we suffer from lattice artifacts due to small N t . Therefore, the 
fixed-scale approach is complemental to the fixed approach. A good feature of the fixed-scale 
approach is that the computational cost for zero-temperature simulations, which is a big burden in 
the fixed approach, can be largely reduced. 

Secondly, the choice of lattice action has important implications for the lattice artifacts. Al- 
though recent large scale simulations of finite temperature/density QCD mostly adopt staggered- 
type lattice quarks [Q, |5|], a rigorous proof is missing for the staggered-type quarks about the ex- 
istence of the continuum limit with desired number of quark flavors. We thus adopt Wilson-type 
quarks with which the continuum limit is guaranteed for any number of flavors. To remove 0(a) 
lattice artifacts, we adopt a clover-improved Wilson quark action |§] with the clover coefficient csw 
non-perturbatively evaluated by the Schrodinger functional method [0]. Note that the Schrodinger 
functional method also requires fine lattices. 

Thirdly, EOS is sensitive to quark masses and to the number of flavors. Therefore, for a 
realistic calculation, we should incorporate the strange quark and try to extrapolate to the physical 
quark mass point. Even an effect of the charm quark on EOS has been recently discussed [j|]. 
With Wilson-type quarks, however, systematic calculation of EOS has been limited to the case of 
two-flavor QCD with rather heavy quarks [Q]. 

We extend the study to the case of 2+1 flavor QCD. The computational cost can be in part 
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reduced by adopting the fixed-scale approach by borrowing zero-temperature configurations from 
large-scale spectrum studies. We adopt configurations of 2+1 flavor QCD with non-perturbatively 
improved Wilson quarks, generated by the CP-PACS+JLQCD Collaboration ^ [l0| ]. The con- 
figurations are available on the ILDG. Using the same coupling parameters, we carry out finite 
temperature simulations varying N t . Details of the simulation parameters are discussed in Sect.|3[ 



2. T-integration method 

In ref . |Q] , we have developed the 'T-integration method" to evaluate the pressure non-perturbatively 
in the fixed-scale approach: Using a thermodynamic relation valid at vanishing chemical potential 

±(P_\e-3p 

dr W r 4 ' K ' 

we obtain 

P f T ,„e-3p 



, dT ^- (2.2) 

T 4 J To T 5 

with p(T ) « 0. Here, the trace anomaly e — 3p is calculated as usual at each temperature. In 
the fixed-scale approach, T is restricted to discrete values due to the discreteness of N t . For the 



integration of (|2.2|), we need to interpolate the data with respect to T. The systematic error from 
the interpolation should be checked. Note that, because the scale is common for all data points in 
the fixed scale approach, T is determined without errors besides the common overall factor I /a. 

3. Lattice setup 

As the zero-temperature configurations, we adopt the results of a 2+1 flavor QCD spectrum 



study with improved Wilson quarks by the CP-PACS+JLQCD Collaboration [jlOp. The QCD action 
S = S g + S q is defined by the RG-improved gauge action S g and the clover-improved Wilson quark 
action S q , 

S g = -A £ c W^ l (x)+ £ erf* 2 ®), (3.1) 
S 9 = £ E«£W, (3-2) 

f=u,d,s x,y 

n 

-S x , y c SW Kf £ o^yF^v, (3.3) 

H>v 

where csw is non-perturbatively determined as a function of /3 [0]. Among the simulation points 
by the CP-PACS+JLQCD Collaboration, we choose j3 = 2.05, K ud = 0.1356 and K s = 0.1351 
which correspond to the smallest lattice spacing a ~ 0.07fm, and the lightest u and d quark masses 
m n jm p ~ 0.63 and m^Jm^ ~ 0.74 in the study. The lattice size is 28 3 x 56 and the statistics is 
about 6000 trajectories. 
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Figure 1: Left: Estimated temperatures of our lattice setup for each N t at j3 = 2.05 and 1.90. Physical scale 
is determined by Sommer scale ro = 0.5fm. Right: The Polyakov loop time history as a status of finite 
temperature Nf = 2+1 QCD simulation with improved Wilson quarks. 



Using the same coupling parameters as the zero-temperature simulation, we are generating 
finite temperature configurations on 32 3 x A/, lattices with N t = 4, 6, ■ • •, 16. Setting the lattice 
scale by ro = 0.5fm, the temperatures for these N t at j3 = 2.05 are shown in the left panel of Fig. [[]. 
The pseudo-critical temperature is expected to be N t ~ 14. Current status of our finite temperature 
simulation is shown in the right panel of Fig.[TJ. Here 1 trajectory is equal to 0.5 molecular dynamics 
step. 

The trace anomaly (e — 3p) /T 4 for our action is given by 

e-3 P tif a W/to\ +a d J^L(4L\ + «**/*\ 1 (3.4) 



T 4 Nf\ da\dp/ sub da\dK ud / mh da\dK s/ mh/ 

||\ = NgN t (-( £ c Q wl*\x)+ £ c lW ^ 2 ( X )) 

P I \ \X,fl>V X,jU,V / 

+ N f^T«f( ^r^a, v F, v (D-%\\ (3.5) 
+c sw /^ £ Tt^g^ v (D-%X\ (3.6) 



dKf 



where (• • -) su fy means that the T = value is subtracted. To evaluate the traces in (3.5) and (3.6) we 



apply a random noise method with complex U(l) random numbers [11]. The number of noise for 
each of the color and spinor indices is 1 . 

4. Beta functions 

The beta functions are required to calculate the trace anomaly (|3T4|). The beta functions are 
obtained through the coupling parameter dependence of zero-temperature observables on a LCP. In 
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Figure 2: The global fit for coupling parameters, Left: j3, Center: K ut i, Left: K s , as the function of m p a. 
Square symbols show coupling parameters in CP-PACS/JLQCD study. The Solid lines show the global fit 
results for each simulation point with corresponding m p /m K and m^/m*. 



our previous report [J|], we have tried to evaluate the beta functions by an inverse matrix method 
It turned out that, although the results for EOS are consistent with an expectation from the 

two-flavor case, errors in the beta functions for the hopping parameters are large. To resolve the 

problem, we adopt an alternative method, the direct fit method [g], to estimate the beta functions. 

We use the results for the hadron spectrum by the CP-PACS+JLQCD Collaboration [10] at three 

jS's, five few's and two jc/s, i.e. totally 30 simulation points at T = 0. 

We fit the coupling parameters, j8, k U( j and K s as a polynomial function of three observables 

am p , m n jm p and m^Jm^ up to the 2nd order, 




2 , r* 



co + ci{am p ) +C2{am p ) +c 3 — 



+ c 4 — +c 5 {am p ) — 
\m p J H \m p 



2 



+c 6 (^)+c 7 (^) + csK)P]+c 9 (^](^). (4.1) 
\m$ J V m / V m <? / \ m p J \ m <t> J 

Figure § shows the results of the global fits ( |Q| ) as functions of m p a. The fits lead to reasonable 
# 2 /dof ~ 1 except for the j3 fit whose ^ 2 /dof is about 5. 

We define LCP by fixing m n jm p and m^Jm^. Therefore, in (4J_), the lattice spacing depen- 
dence on a LCP appears through the terms containing am p . The beta functions are thus calculated 
from the coefficients c\, ci, C5 and c§. From the fits we obtain preliminary values for the beta func- 
tions: a^- = -0.334(4), J^- =0.00289(6) and a^- = 0.00203(5) at our simulation point. 

aa aa aa 

Here the errors are statistical only. Estimation of systematic errors is left for future investigations. 



5. Equation of state 

The left panel of Fig. |3| shows the results for the trace anomaly, together with its decomposition 
into /? and K derivative parts in (^4|). The lines are spline interpolations. We find that there is a 
big cancellation between the j3 and K derivative parts and the resulting peak height of the trace 
anomaly is about 7. This relatively low peak height obtained at N t ~ 14 is roughly consistent with 
recent results from highly improved staggered quarks obtained on N t = 6-12 lattices using the 
conventional fixed N t approach [Q|5]]. 
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Figure 3: Left: The trace anomaly and its j3 and jc's derivative contributions are shown. The latter contribu- 
tion is a minus quantity. The lines are drawn by a spline interpolation. Right: The trace anomaly and EOS, 
e/T 4 and 3 p/T 4 are shown. The lines are drawn by a spline interpolation. 



Carrying out the T-integration CL2) using a trapezoidal interpolation of the trace anomaly, we 
obtain the pressure p/T 4 shown in the right panel of Fig. |3[ Here, we have chosen the starting 
point of the integration to be at N t = 16 where the trace anomaly vanishes within the statistical 
error. The energy density e/T 4 is calculated by p/T 4 and (e — 3p)/T 4 . The lines in the figure are 
spline interpolations. 

The overall large errors in p/T 4 and e/T 4 are propagated from the large errors in (e — 3p) /T 4 
at low temperatures in the numerical integration. From the left panel of Fig. ^j, we note that the 
large error in (e — 3p) /T 4 at T ~ 200 MeV is due to the /3 derivative part in (JO]). We are currently 
trying to increase the statistics on N t = 12-16 lattices. 

At the same time, we are starting simulations at jS = 1.9 to study the scaling. Corresponding 
temperatures at j3 = 1.9 are shown in the left panel of Fig. p]. Smaller N t in the low temperature 
region will make the statistical error problem less severe. Another objective to study at /3 = 1.9 is 
to extend the investigation of EOS to lighter quark masses [12]. 



6. Summary 

We presented the status of our project to calculate EOS for 2+1 flavor QCD with improved 
Wilson quarks. Previous studies with Wilson-type quarks were limited to the case of two-flavor 
QCD with not quite light quarks [||]. Adopting the fixed-scale approach, we succeed to calculate 
the first EOS for 2+1 flavor QCD with Wilson-type quarks on a fine lattice. Although the light 
quark masses are heavier than the physical values yet, our EOS looks roughly consistent with 
recent results with highly improved staggered quarks [|], g]. On the other hand, we note that a 
large cancellation in the zero-temperature subtraction of the gauge action in the low temperature 
region leads to large statistical errors in the final EOS. We need a large statistics there. However, 
because the problem is limited on several lattices, we think that the overall computational cost is 
much smaller than that required with the conventional fixed N t approach, and tractable with current 



6 



EOS in 2+1 flavor QCD with improved Wilson quarks 



T. Umeda 



computer powers. We are now ready to start a study with much lighter quarks, adopting the on-the- 



physical -point configurations by the PACS-CS Collaboration [ |12| ] 
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